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Abstract 

Starting from Martin, Roberts and Stirling fit for unpolarized deep 
inelastic structure functions and using the newest experimental data 
on spin asymmetries we get a fit which provides polarized quark dis- 
tributions. We analyze the behaviour of such functions near x equal 
to 1. The first moments of these distributions are also discussed. Our 
fit prefers combination of proton and neutron data versus proton- 
deuteron one. 
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Recently the interest in the spin structure of the nucleon has been re- 
newed due to new data taken in experiments at CERN [[l], |2[ and SLAC 
[^]. The polarized deep inelastic asymmetries for deuteron and proton (Spin 
Muon Collaboration at CERN) as well as for neutron (at 3 He target in E 
142 experiment at SLAC) were precisely measured also in small x region. 
Together with an old SLAC |4[] and EMC || data for proton one has con- 
siderable amount of information which can be used to study nucleon spin 
structure and particularly quark parton distributions. 

The unpolarized quark distributions in the nucleon are known due to 
several fits |6|, |7| to the experimental data. The one of the most recent 
fit is given by Martin, Roberts and Stirling (MRS) || who used existing 
experimental points in order to determine quark and gluon distributions. 

In this paper we would like to determine polarized quark parton distri- 
butions starting from the MRS fit and using all existing data for proton, 
neutron and deuteron spin asymmetries. We shall concentrate not only on 
a controversial low x region but also on the behaviour of polarized parton 
distributions at x — > 1. The phenomeno logical analysis of CERN and SLAC 
data gives different results for AS (quark spin content of the proton) and 
As (strange quark polarization in the proton). Our results, as one will see 
later, prefers rather the numbers gotten from an analysis based on neutron 
E142 data. 

Let us start with the formulas for unpolarized quark parton distributions 
(at Q 2 = 4 GeV 2 ) given by Martin, Roberts and Stirling. First of all we shall 
consider their fit called D'_ with very singular behaviour of " sea" distribution 
at small x values (which, however agrees with the results from HERA ||). 
We have for the valence quarks distributions: 

u v (x) + d v (x) = 1.422aT - 58 (l - xf 92 (l + 2.59y/x~ + 4.21a;), 

d v (x) = 0.074x~ a76 (l-x) 4 - 67 (l + 28.7 v / x + 8.58x), (1) 

and for the "sea" ones: 



2u(x) = 0AS(x)-5(x), 

2d{x) = 0AS(x) + 5{x), (2) 
2s(x) = 0.2S{x), 
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where 



S(x) = 0.083ar L5 (l - x) 7A (l + 8.57v^ + 15.8s), (3) 

and 

5(x) = 0.164x~ a58 (l-x) 7 - 4 . (4) 

The most important features of this fit are that the valence "up" quark 
distribution dominates at x — > 1 and that the "sea" is not SU(2) symmetric 
(having Lipatov type behaviour || at small x). Because the unpolarized 
parton distributions are the sum of spin up and spin down distributions 
whereas the polarized ones are the difference of these functions our idea is 
just to split the numerical constants in formulas (1,3,4) in two parts in such 
a manner that we get positive defined distributions (which is not so easy to 
achieve). Our expressions for Aq{x) = q + (x) — q~(x) (q(x) = q + (x) + q~(x)) 
are: 



Au v (x) + Ad v (x) = x~°- 58 (l - x) 3 - 92 (a 1 + a 2 \fx + a 3 x), 
Ad v (x) = x- a76 (l-x) 4 - 67 (6 1 + 6 2V / S + M, 

(5) 

AS(x) = cx- a58 (l-x) 7 - 4 , 
A5(x) = dx-°- 58 (l-x) 7A . 

We want to stress that we are using very simple way of parameterization 
of q + and q~ without introducing any additional powers of x or (1 — x), so 
only the coefficients for polarized quark distributions are fitted. In order to 
get the integrable function AS(x) we divide two first coefficients in eq.(3) 
into equal parts. 

We fit our formulas (with eight parameters) to the experimental data on 
spin asymmetries, given by: 



A\{x) 



Al{x) 



AAu v (x) + Ad v {x) + 2.2AS{x) - 3A<5(x) 
Au v {x) + djx) + 2.2S(x) - 35 (x) 

A^(ar) + AAd v (x) + 2.2AS(x) + 3A5(x) 
ujx) + 4djx) + 2.2S{x) + 35{x) 



(6) 
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In this paper we assume that the spin asymmetries do not depend on Q 2 
what is suggested by the experimental data |], |3| and phenomenological anal- 
ysis [ IT5 1 . The unknown parameters in eq.(5) are determined by making best 



fit to the measured spin asymmetries for proton (SLAC-Yale, EMC, SMC) 
and neutron (E142). We get the following values (taking care of positivity 
for quark distributions): 

a x = 0.874, a 2 = 5.023, a 3 = 12.73, 

6i = 0.074, b 2 = 0.884, b 3 = 0.649, (7) 

c = -0.556, d = -0.004. 

It is interesting to note that the fit shows no significant SU(2) symmetry 
breaking for "sea" polarization (d coefficient is close to zero). The x 2 is 20.1 
for 34 degrees of freedom for such a fit. In figures la, lb and 2 the comparison 
of our fit with the experimental points for proton and neutron asymmetries 
is given. The prediction for the deuteron case (which adds 7.2 to x 2 f° r H 
degrees of freedom) with Af given by the following expression: 

d _ 5Au v (x)+4AAS(x) 3 
A{X) ~ 5u v (x)+4AS(x) (1 " 2 Pd) ' (8) 

(Pd is D-state probability equal to 5.8%) is presented in figure 3. We see 
that in this case the calculated curve at small x lies above the experimental 
points. This is because we get positive Af for all values of Bjorken variable, 
whereas the data are negative in small x region. Of course we can make a 
fit with inclusion of SMC deuteron data. Than we get x 2 = 26.6 (for 45 
degrees of freedom) that is not much less than for proton+neutron case and 
predicted deuteron asymmetry (x 2 = 20.1 + 7.2 = 27.3). Making the fit 
to proton+deuteron asymmetries we get x 2 = 22.2 (37 degrees of freedom) 
which become 41.3 when one adds the result for A™ (x 2 = 19.1 for eight 
points). Hence, it is possible to get a satisfactory deuteron asymmetry using 
proton+neutron data and is not when one tries to get the neutron asymmetry 
from proton+deuteron data. This is our x 2 argument why in our fit we use 
SLAC neutron data omitting CERN deuteron ones. As we will see below 
such preference is justified also when one analyzes an integrated quantities. 

We can use our fit to determine the first moments of parton distributions. 
We get e.g. for I = ft g x {x, Q 2 ) dx at Q 2 = 4 GeV 2 : 
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14 1 
7 n = — Au + — Ad + — As = -0.027. 
18 18 18 



(9) 



Other combinations (singlet and octet ones) of quark polarizations are: 



a = AS = Au + Ad + As = 0.50, 

a 3 = Am - Ad = 1.23, (10) 

a 8 = Au + Ad - 2As = 0.71, 

whereas As = —0.07. 

We would like to stress that these results are obtained from the fit without 
any constraints (e.g. for integrated quantities), which is not a case in other 
fits (see e.g. ref.[0). We also do not use any information from neutron and 
hyperon /3-decays. The value of a 3 should be compared to qa = 1.257 ±0.003 
[0 lowered by several percent due to QCD corrections (see e.g. ref.0]), i.e. 
to 03 ~ 1.1. The agreement is not perfect, although satisfactory. Also the 
value of I p , when compared with the last published number: I p = 0.136 ± 
0.011 ± 0.011 [TJ], is too high. This happens because the structure function 
Fx (x) given by the MRS fit lies few percent higher (in considered region) than 
the one used in phenomenological analysis by EMC and SMC. Because the 
ratio of unpolarized to polarized structure functions, i.e. the asymmetry is 
fitted very well the result for g\{x) is also overshooted. The second difference 
arises because the experimental groups use " Regge" type extrapolation in the 
x —>■ region. For example, the latest SMC value for I p is gotten assuming 
constant value for gi(x) in unmeasured small x region. When we modify our 
distributions in such a manner (only for x between and 0.003) we get: 

I p = 0.168, J n = -0.038, AS = 0.43, , . 

a 3 = 1.23, a 8 = 0.62, As = -0.06, ^ ' 

the values which are not substantially different from those presented in 
eqs.(9,10). The second reason for using this modification (for small x) is that 
MRS distribution for u v (x) is not positive defined for tiny x's (x ~ 10~ 7 ) and 
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hence such blunder is also present in our polarized distributions for u+(x) 
and u~(x). 

If we fit the polarized parton distributions to all measured asymmetries 
(on proton, neutron and deuteron targets) we get e.g. a 3 = 1.50 which is 
approximately 40% to big. If one uses D Q fit of Martin, Roberts and Stirling 
instead of D_ (D is not so divergent at x ~ 0) we get the similar values 
for integrals I p and I n , whereas 03 = 1.34. Also the x 2 is worse in this case. 
Hence, we prefer our distributions fitted to asymmetries measured on proton 
(SLAC-Yale, EMC, SMC) and neutron (SLAC E142) targets and having its 
roots in the D_ fit. 

Now, we would like to make some comments about the x — > 1 behaviour 
of valence quark distributions. Looking at the data points for proton spin 
asymmetry the value close to 1 at x ~ 1 is preferred, whereas for neutron 
and deuteron case values close to seem to be natural (such observation is 
fragile due to the big experimental errors in this x region). In our approach 
we can give predictions for the behaviour of polarized quark distributions 
and spin asymmetries in the x — > 1 limit. Let us assume that valence u and 
d quark polarized distributions behave at x — > 1 as: 



v>v — > a±(l - x) p + . . . , 

d ± ^ b ± (i-xy + ..., (12) 

where p is the smallest power of (1 — x) term. Than, we have in such a limit: 

F? a + -a_ +4(6+ - &_) 
Ff ~* 4(a+ + a_) + b+ + 6_ 

and: 



(13) 



A\{x) 



4(a+ 



4(a + + a_) + 6+ + 6_ 



1V 1 a + + a_ +4(6+ + 6_) v 1 

In the case of SU(6) symmetry the flavour-spin part of nucleon wave 
function gives: a+ : a_ : 6+ : 6_ = 5 : 1 : 1 : 2. Hence, one gets well known 
results: F"/Ff — > 1/4, whereas A\ — > 5/9 and A™ — > 0. The counting rules 
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for parton distributions at x — > 1 (see ref . ||14|| ) yield: a^/a + = 6-/6+ = 



so one has A\,A\ — > 1. Authors of ref. [14| assume in addition: a + /b + = 5 
(as in the SU(6) symmetric case), hence they get F™ / Ff — > 3/7. In the 
MRS fit one has 6±/o± = which leads to Fp/Ff -> 1/4. For the spin 
asymmetries one gets A\ = A\ (a+ — a_)/(a+ + a_) < 1. In our case 
we have A\ = A\ = A\ = 0.77. The different limits (at ac 1) for A? and 
A™ one can get assuming that coefficients 6 are not negligible in comparison 
to a's (see eq.(14)). But then it is impossible to get the value for F™ / Ff 
suggested by the experimental data, namely 1/4 (see the figures in ref. [|iT| ) . 

We have not considered gluon distributions in our fit (there is no need 
to take it into account dealing with functions at fixed Q 2 ). If, however we 
include explicit gluon terms into the asymmetries (in the way proposed in 
ref. ||15|| ) we get a fit which is worse than the considered one. 

Starting from the MRS fit B to the unpolarized deep inelastic scattering 
data we made a fit to proton and neutron spin asymmetries in order to obtain 
polarized quark parton distributions. We got Am = 0.91 (Au v = 1.04), Ad = 
—0.33 (Ad v = —0.18) and As = —0.07 for integrated quantities. At x — > 1 
the asymmetries for nucleons point towards the value equal to 0.77. With 
the improved behaviour at x — > for the unpolarized parton distributions 
and consistent data for spin asymmetries (mainly for deuteron) our method 
of determination of quark polarized distributions looks promising. 
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Figure captions 

Figure la The comparison of spin asymmetry on protons (data points from SLAC- 
Yale, EMC and SMC experiments) with the curve gotten from our fit 
(eqs.(5,6,7)). 

Figure lb The same as in figure la but with x in logarithmic scale. 

Figure 2 The comparison of spin asymmetry on neutrons (SLAC E142 data) 
with the curve gotten from our fit (eqs.(5,6,7)). 

Figure 3 Our prediction for deuteron asymmetry (NMC data). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409204vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9409204vl 



Figure lb 




This figure "figl-3.png" is available in "png" format from: 
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Figure 2 
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